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ABSTRACT 

The large scatters of luminosity relations of gamma-ray bursts (GRBs) have been one of the 
most important reasons that prevent the extensive applications of GRBs in cosmology. Many 
efforts have been made to seek tight luminosity relations. With the latest sample of 116 GRBs 
with measured redshift and spectral parameters, we investigate 6 two-dimensional (2D) correla- 
tions and 14 derived three-dimensional (3D) correlations of GRBs to explore the possibility of 
decreasing the intrinsic scatters of the luminosity relations of GRBs. We find the 3D correlation 
of -Bpcak-TRT-i to be evidently tighter (at the 2a confidence level) than its corresponding 2D 
correlations, i.e., the Ep^-L and trt-L correlations. In addition, the coefficients before the 
logarithms of i? pea k and trt in the E pea ^-T^n-L correlation are almost exact opposites of each 
other. Inputting this situation as a prior reduces the relation to L oc (£ ; poa k/' r RT) ' 842±0 ' 064 j 
where E peak and r RT denote the peak energy and minimum rise time in the GRB rest frame. We 
discuss how our findings can be interpreted/understood in the framework of the definition of the 
luminosity (energy released in units of time). Our argument about the connection between the 
luminosity relations of GRBs and the definition of the luminosity provides a clear direction for 
exploring tighter luminosity relations of GRBs in the future. 

Subject headings: gamma-ray burst: general 
1. Introduction 



Gamma-ray bursts (GRBs) are the most lumi- 
nous astrophysical events observed so far. The 
high luminosities make them observable at very 
high redshifts. The measured highest redshift 
for GRBs has exceeded 8 (e.g ., the GRB 090423 
with a redshift of z w 8.2 dTanvir et all 120091: 
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Salvaterra et~al1l2009h and the G RB 090429B with 
a pho tometric redshift of z ~ 9.4 (jCucchiara et al 
20111 )). Many efforts have been made to apply 



them to cosmology (see, for example. iDai et al.l 
2004: Ghirlanda et alJl2004ri iFirmani et al.ll200£ 
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is to find some luminosity relation that relates 
the luminosity (e.g., the isotropic peak lumi- 
nosity L) or energy (e.g., the isotropic energy 
£^,iso or the collimation-corrected energy E 7 ) of 
GRBs to their measurable properties (the lumi- 
nosity indicators). The tighter the relation is, 
the more accurate and reliable information about 
the universe we can get from it. Many luminos- 
ity relations have been discovered in the liter- 
ature , e.g., the relations of n aa -L dNorris et ah 
2000), V-L (|Fenimore fc Ramirez-Ruij [2000: 
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Reichart et al, l200lh. ■E n eak~-fi'-Y.isn (lAmati et al 
20021). E^yE-t (iGhirlanda , et al.l 2004ah . E^~- 
L dSchaefeijl200ilWei fc Gaol2003tlYonetoku et al 
20041 ). and trt-L (|Schaeferll2007l) . etc. Here, the 
spectral lag Ti ag is the time shift between the hard 
and soft light curves. The variability V is a quan- 
titative measurement of the spikiness of the light 
curve, which can be obtained by calculating the 
normalized variance of the observed light curve 
around the smoothed light curve. There exist 
several definitions of V , depending mainly on the 
smoothing time intervals upon which the reference 
curve is built and the normalization. -Epoak is the 
photon energy at which the vF v spectrum peaks. 
The minimum rise time trt of a GRB light curve 
is the shortest time over which the light curve rises 
by half the peak flux of the pulse. In addition, cor- 
relations were also found between the transition 
times of the X-ray light curve from exponential 
to power la w and the X-ray luminosit i es at the 
trans itions (|Dainotti et all l2008l l2010t iQi fc Lu 
2010h . 

The most crucial problem with GRBs for cos- 
mological application is that their luminosity re- 
lations are usually quite scattered, which prevents 
them from being used as good standard candles 
as Type la supernovae. Therefore, one of the 
most important directions of cosmological study 
for GRBs is to seek GRB luminosity relations 
with sufficiently small intrinsic scatters. Many 
attempts have been made in literature, among 
which, a remarkable approach is to explore pos- 
sible hidden parameters in known correlations in 
ord er to reduce the intri nsic scatters. For exam- 
ple, Firmani et al. ( 20061 ) claimed that a temporal 
parameter of the prompt emission, To. 45, could re- 
duce the scatter of the correlation of E vca ^-L to a 
negligible value. It was later found that the new 
proposed r elation does not appear to be as tight as 
it seemed ( Rossi et al. _ 2008 ;_ Collazzi fc Schaefer 
20081 ). iTsutsui et al.l (|2009l boioh suggested that 



introducing the parameter Tj, = E~ 1 _i SO /L to the 
Epeak-L relation could substantially reduce the 
intrinsic scatters of the correla t ion o f E pea k-L 



and E pea k~E 1: 

iso • 

IXu fc Huanj (|2012l ) reported 



that a significantly tighter correlation can be ob- 
tained by adding E~, 



ered in Dainotti et al 



to the correlation discov- 



(|2008l . l2010D . However, they 
did not perform the normalization (see the follow- 
ing section of this paper for details about nor- 



malization) before comparing the intrinsic scat- 
ters. Taking into account the normalization, the 
new relation is in fact looser instead of tighter in 
th e sense o f cosm ological distance measurement. 
In Yu et al. I (l2009l) . more general analysis was car- 
ried out, investigating 5 two-dimensional (2D) cor- 
relations and 10 derived three-dimensional (3D) 
correlations to explore the possibility of reducing 
the intrinsic scatters. 



Since we have more G RBs at hand now ([Xiao fc Schaefer 
20091: IWang et al.ll201l[ ). it is beneficial to update 
the comparison between 2D and 3D luminosity 
relations of GRBs and to check whether the con- 
clusions drawn from an earlier GRB sample still 
hold. This is the main aim of this paper. We 
included more GRB lu minosity r e lation s in our 



Six 2D 



analysis than there are in lYu et al. 
correlations and 14 derived 3D correlations were 
investig ated. We basic ally followed the method 
used in lYu et aL ( 20091 ) for the analysis. To be 
complete, we briefly describe the method in the 
following section. The results and discussion are 
presented in Section 3. A summary is given in the 
last section. 

2. Methodology 

For the 2D luminosity relations, we consider the 
following correlations: 
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The first five of the above correlations a re exactly 



the lu minosity relations considered in lYu et al 
(|2009f) . With the latest compilation of GRBs, the 
V-L relation has becom e quite scattered (jXiao fc Schaefer 
20091: IWang et al. 1 120111) . which makes the param- 
eter V a very poor luminosity indicator. Despite 
this, we still include this relation in our analysis in 
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order to investigate whether V is helpful in reduc- 
ing the intrinsic scatter when constructing 3D lu- 
minosity relations together with other luminosity 
indicators. In this paper, we added in our analysis 
the correlation of Eq. ©, for which more GRBs 
can be utilized, compared with the similar corre- 
lation of Eq. dU, due to its independence from the 
GRB jet opening angle. In the correlations, the 
isotropic peak luminosity L, the isotropic energy 
E 7l i SO , and the collimation-corrected energy E^ 
are derived from the observables of the bolometric 
peak flux Pboio, the bolometric fluence SWo, and 
the beaming factor -Fbeam through 



,( 4 ) 



L = 47rd L Pboio, 



E 



7, iso 

7 



= 47rd i 5boio(l + z) 

'7,iso^ 



(7) 

(8) 
(9) 



where is the luminosity distance, which de- 
pends on the cosmological model and is inversely 
proportional to the value of the Hubble parameter 
of today. As in Yu et al. ( 20091 ). we adopt the flat 
ACDM model with J7 m = 0.27 and replace dh with 
i £ = sa.d t x 1 cm in the calculation, so that the 
dependence on the Hubble constant is absorbed 
into the intercepts aj of the linear luminosity re- 
lations. 

The method used to extend the 2D luminosity 
relations to 3D ones is also basi cally the same as 
that adopted in lYu et al.l (2009), except we have 
more 3D luminosity relations here due to the ad- 
dition of the correlation of Eq. ([6]). First, for later 
convenience when denoting specific correlations, 
we summarize the equations for the 2D luminos- 
ity relations above as 
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The coefficients c are given two superscripts to in- 
corporate the 3D correlations introduced below. 
We denote a luminosity relation by the superscript 
pair of the corresponding c. We let i < j for all 
the 3D luminosity relations (i, j) to avoid duplica- 
tion, and we classify them into the following three 
groups for convenience: 

1. Correlations between the luminosity (the 
isotropic peak luminosity L ) and two lumi- 
nosity indicators. For (i, j) with both i and 
j in (1, 2, 3, 5), the 3D luminosity relations 
are 

y^ =cf 3) +cf 3) x^ +c^ 3) x^. (19) 

2. Correlations between energy (the isotropic 
energy Ej^ so or the collimation-corrected en- 
ergy E-y) and two luminosity indicators. For 
(i, j) = (1, 4), (2, 4), or (4, 5), the 3D lumi- 
nosity relations are 

y (A) = c {i,i) + c p) x W + cpVi), (20) 

and for (i, j) = (1, 6), (2, 6), or (5, 6) they 
are 

(6) _ J.i,3) , , r (iJ) v U) 



y K ' = cq 



K {' J >x^> +c%' J1 xU>. (21) 



3. Correlations between the luminosity, the en- 
ergy i? 7i iso or E^, and the peak energy -E pe ak- 
For (i, j) = (3, 4) or (3, 6), the 3D luminos- 
ity relations are 

y«= c P+cf' J V 4 >+cpy^. (22) 

A luminosity relation can be multiplied by an 
equal constant on both sides of the equation with- 
out actually changing the correlation itself. How- 
ever, the multiplication would change the intrin- 
sic scatter of the correlation by a factor of the 
absolute value of the constant multiplied on the 
equation. To compare the intrinsic scatters of 
the luminosity relations, we divide both sides of 
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Eq. ((221) b y a factor of 1 - c 



so that log(di) 



terms have the same coefhcient in all the luminos- 
ity relations discussed here. We normalize the lu- 
minosity relations in this way because, in addition 
to being helpful in understanding the GRBs them- 
selves, the luminosity relations are mainly aimed 
at distance measurements. 

To explore possible hidden parameters in the 
2D luminosity relations, we compare the intrin- 
sic scatters of 3D luminosity relations with those 
of corresponding 2D ones. The principle is that 
the two luminosity relations compared with each 
other should have two parameters in common, so 
that we can conclude whether the intrinsic scat- 
ter is reduced by introducing the third parame- 
ter. Following this principle, for the 3D luminos- 
ity relations (i, j) in the first and third classes, 
i.e., those of Eq. (fl9l) and (f22j) . they are compared 
with 2D luminosity relations (i, i) and (j, j), and 
for the 3D luminosity relations in the second class, 
i.e., those of Eq. (f20f and (l2~T1) . they are compared 
with 2D luminosity relations (4, 4) and (6, 6) re- 
spectively. Only when the intrinsic scatter of a 
3D luminosity relation is smaller than that of all 
its corresponding 2D luminosity relation(s) can we 
say that the intrinsic scatter is reduced. 

In the fit of the lumin osity relations, we used 
the techniques presented in lD'Agostinil(|2005h . Us- 
ing the correlation of Eq. (|19|) as an example, the 
joint likelihood function for the coefficients c and 
the intrinsic scatter <7; n t is given as 
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(23) 



where k runs over GRBs with corresponding quan- 
tities available. The joint likelihood functions are 
similar for other 3D luminosity relations. For the 
2D luminosity relations, the joint likelihood func- 
tions can be obtained by setting c 2 = 0. 

For the GRB data, we used the compilation 
which includes 116 GRBs. 



Wang et al 



(120111 ). 

For this compilation, th e data of the luminosit y 
indicators are taken from lXiao fc Schaeferl (|2009f h 
The size of this GRB sample is almost twice as 



large as the compilation in I Schaeferl (|2007l) . which 
was used in the previous s tudy to compar e 2D and 



3D luminosity relations (|Yu et al. 2009). In ad 



dition, significant improvements have been made 
in the calculation of the luminosity indicators for 
the updated compilation. When considering error 
propagation from a quantity, say £ with error <jf, 
to its logarithm, we set [log(£+cr^)+log(£— a^)]/2 
and [log(£ + <r^) — log(£ — a^)]/2 as the center 
value and the error of the logarithm correspond- 
ingly. This requires £ > (the quantities we are 
interested in here are all positive) . Due to the lim- 
itation of the data, for a given luminosity relation 
(i, j), not all the GRBs have all of the needed ob- 
servational quantities available and satisfy £ > cr^T 
at the same time. By set (i, j) we denote the max- 
imum GRB set that can be used in the luminosity 
relation (i, j). The numbers of GRBs of different 
sets are presented in Table [TJ 

3. Results and discussion 

We summarize our results for the fits to the 
luminosity relations and comparisons between 2D 
and 3D correlations in Table [T] We list our find- 
ings from the table as follows: 



1. 



For the c orrelation s com mon with those 
studied in lYu et all (|2009h . including both 
2D and 3D ones, almost all the intrinsic scat- 
ters have increased with the updated GRB 
sample. For the 2D ones, only the tkt-L re- 
lation retains a comparable intrinsic scatter. 

2. For the cases of (1, 3) and (1, 5), the in- 
trinsic scatters of the 3D correlations are re- 
duced at the lcr confidence level compared 
with their corresponding 2D correlations, 
and for the case of (3, 5), the intrinsic scat- 
ter is reduced at the 2a confidence level. No 
other statistically significant reduction of the 
intrinsic scatters is found. 

3. For the case of (3, 5), the values of the coef- 
ficients Ci and C2 are almost exact opposites 
of each other. For the case of (1, 3), the val- 
ues of ci and c 2 are also rough opposites of 
each other. 

Item 1 reminds us once again of the complexity 
of GRBs. Correlations should be tested against a 
large enough sample to determine how good they 
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Table 1: Fit of 2D and 3D luminosity relations. The first row and the first column of the table are the indices denoting the luminosity 
relations and the corresponding luminosity indicators. In every cell of the table, the first row is the luminosity and/or energy involved in the 
relations which, together with the luminosity indicators, tell what the correlations are about. The second row of every cell is the number of 
GRBs of set (i, j), the vector below enclosed by parentheses is the vector of c for the luminosity relation (i, j), and what follows next is the 
intrinsic scatter. For 3D luminosity relations, the reduction of the intrinsic scatters compared with corresponding 2D luminosity relations 
is presented in the brackets. The statistics in the table are for the median values and the errors with the la (68.3%) confidence level. 



are. The intrinsic scatter derived from a small 
sample may be significantly affected by the selec- 
tion effect of the sample itself. Item 2 is basically 
consis tent with the previous studies of lYu etal 



(l2009h . in which the intrinsic scatters are reduced 
at lcr and 2a confidence levels, respectively, for 
the cases of (1, 3) and (3, 5), i.e., the Ti ag -E pca k- 
L and E poli k-Tm-L correlations. With the new 
GRB sample adopted here, the intrinsic scatter is 
also reduced at the la confidence level for the case 
of (1, 5), i.e., the ti^-t^-L correlation, the edges 
of the two la confidence intervals of the magnitude 
of reduction of the intrinsic scatter for (1, 5) are 
both very close to zero, which reduces its reliabil- 
ity. The correlation (3, 5) turns out to be the only 
3D correlation that appears very robust in the test 
of the reduction of the intrinsic scatter. 

The correlation (3, 5) is among -Epeak, Tit, 
and L. It is a relation among an energy scale, 
a timescale, and the luminosity. This naturally 
makes us think of the definition of the luminos- 
ity (energy released in units of time). As can be 
seen from the large scatters of the 2D luminos- 
ity relations, it is most likely that, due to their 
complexity, the luminosities of GRBs cannot be 
well determined from only one quantity. How- 
ever, according to the definition of the luminos- 
ity, we should in principle be able to calculate the 
luminosity from a characteristic energy scale and 
a characteristic timescale of GRBs. If we could 
find two measurable quantities (e.g., an energy 
scale and a timescale) that are strongly correlated 
with the assumed characteristic energy scale and 
timescale, then we would find a 3D luminosity re- 
lation between the luminosity and the two mea- 
surable quantities. The quality of the relation ob- 
viously depends on the strength of the correlations 
between the two quantities and the assumed char- 
acteristic energy scale and timescale, and is also 
related to the degree of match between the two 
quantities in the sense of the match between en- 
ergy scale and timescale in the definition of the lu- 
minosity (i.e., the energy scale and timescale cho- 
sen to calculate the luminosity should correspond 
to each other). The intrinsic scatter of such a 3D 
luminosity relation is expected to be significantly 
reduced compared with its corresponding 2D ones. 
The situation described in item 2 suggests .Epcak 
and trt as one such pair of measurable quantities. 



-2 

-2.5 




-1.5 -1 -0.5 0.5 1 

i pW(i+») / mr(i+«r 1 l 

& 300 keV ' 0.1 a 



1.5 



Fig. 1. — Best fit of the GRBs to the luminos- 
ity relation of Eq. (j2~4")l and the corresponding 2a 
confidence region. In the calculation of the L, the 
dependence of the luminosity distance on the Hub- 
ble constant is absorbed into the intercept of the 
relation. 



What makes the correlation (3, 5) more inter- 
esting is the situation described in item 3. Follow- 
ing the discussion above, the almost exact opposi- 
tion of ci and c-i suggests a considerably high de- 
gree of match between £p e ak and trt considering 
that, in the definition of the luminosity, the indices 
for the energy scale and the timescale are oppo- 
site. In comparison, the degree of match between 
-Epeak and Tiag is slightly weaker. The situation of 
the almost exact opposite for the correlation (3, 5 ) 
does not show up in the results of Yu et a l. (2009). 
However, remember that, for the GRB data used 
here, not only is the sample size much larger, but 
significant improvements also have been made in 
the calculation of the luminosity indicators, so we 
have reasons to treat the situation as more than 
just a coincidence. If we input the opposition of 
ci and C2 as a prior for the correlation, then the 
relation is reduced to a 2D one, that is, 
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A fit to this relation gives the result of a = 
-3.742±°;°^, b = 0.8421°;°^ and the intrinsic 
scatter a mt = 0-346to.o32- We present the relation 
in Figure [TJ Using E pcak and r RT to denote the 
peak energy and minimum rise time in the GRB 
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rest frame, the relation can be expressed as 

. 0.842±0.064 
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Following the definition of the luminosity, if appro- 
priate luminosity is used on the left-hand side of 
Eq. (HU or Eq. (g5]), the coefficient b in Eq. 
or the index in Eq. (1231) should be 1. However, 
in practice, there are various definitions of lu- 
minosities, e.g., peak or averaged luminosity and 
isotropic or collimation-corrected luminosity. It is 
unclear which luminosity is attached to the pair of 
energy scale and timescale used on the right-hand 
side of equations, i.e., -Speak and trt- What we 
adopt in the relations is the isotropic peak lumi- 
nosity L. Thus, the coefficient b in Eq. (f2"4"|) or the 
index in Eq. (1251) may be interpreted as the cor- 
rection factor between the luminosity attached to 
the pair of Sp Ca k and trt and the isotropic peak 
luminosity. 

For a clear and complete physical interpretation 
of the £? P eak-TRT-£ correlation, it should be ex- 
plained how Speak and trt are correlated with the 
assumed energy scale and timescale. This is an is- 
sue that needs further detailed study, which is out 
of the scope of this short paper. Here, we present 
some general speculation on this. One may note 
that Speak is not a radiated energy, but the pho- 
ton energy at which the vF v spectrum peaks. The 
calculation of the luminosity needs the radiated 
energy. S pea k may relate itself with a radiated en- 
ergy through a physical process that releases pho- 
tons with the frequency at which vF v spectrum 
peaks, or it may just be a dominant parameter 
in the spectrum that can be used as a represen- 
tation of the bolometric energy released by the 
GRB, as can be seen from E pea ^-E 7 ilso and S P e a k _ 
S 7 correlations (S 7! i so and S 7 cannot replace the 
role of S p c a k here, mostly because their calcula- 
tion depends on the luminosity distance and the 
intrinsic scatters are enlarged in the procedure of 
normalization, while S p0 ak is directly measured 
from the observation and is independent of the 
cosmological model). We think that trt stands 
out in our investigation most likely because trt 
is the minimum rise time. Though requiring it to 
be minimum sounds like a very simple operation 
in the data processing, it could be a very impor- 
tant step for extracting the clean physics from the 



complex GRBs, thus eliminating intrinsic scatters 
from many irrelevant physical processes. 

We would also like to mention the 3D cor- 
relation of (3, 6), i.e., the S P eak _ S 7i i so -S corre- 
lation and the S D -T^-S D correlation referenced 



Tsutsui et all (|2009l . l2010h . Using the notation 



adopted here, the latter is given by 

logi = A + Slog[S P eak(l + z)\ + C\ogT L , (26) 

where Tl — S 7 .i SO /S. Substituting Tl into the 
equation, the correlation can be rewritten as 

(1+C) logS = ^+Sl0g[S p eak(l+2)]+Cl0gS 7 . iso . 

(27) 

Comparing it with Eq. (|2"2"1) . we can see that this 
relation is in fact the same as the 3D luminosity 
relation (3, 6), and it has been already normal- 
ized. However, the conclusion we draw here for 



the correlati o n conflicts with that inlTsutsui et aL 
(|2009ll2010l) . iTsutsui et all (j2009|, |2010|) claimed 



that the intrinsic scatter of the correlation E p -Tl~ 
L p is smaller than that of both the S pea k _ S and 
S P eak~S 7! i so correlations. Our result shows that 
the intrinsic scatter of the correlation of (3, 6) is 
not reduced compared with its corresponding 2D 
correlations. In fact, its intrinsic scatter even ap- 
pears a little larger. The advantages of our analy- 
sis here are that we used more GRBs, which helped 
us reduce the selection effect from the sample it- 
self, and that we adopted a sophisticated statisti- 
cal method, in which not only the intrinsic scatter 
but also its error can be derived, so that we can 
judge the statistical significance of the reduction 
of the intrinsic scatters. 

4. Summary 

As a further step on in seeking tight luminos- 
ity relations of GRBs, we fitted the latest data of 
GRBs with measured redshift and spectral param- 
eters to 6 2D correlations and 14 derived 3D cor- 
relations and compared their intrinsic scatters to 
explore possible hidden parameters in the 2D cor- 
relations. Compared with the analysis of an early 
sample of GRBs, the intrinsic scatters of most of 
the luminosity relations have increased, which re- 
minds us of the complexity of GRBs. Correlations 
should be tested against a large enough sample to 
reduce the possible selection effects from the sam- 
ple itself. Our result shows that the S pC ak-7RT-S 
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correlation appears to be significantly tighter (at 
the 2a confidence level) than its corresponding 2D 
correlations, i.e., the E pe3 ^-L and trt~L correla- 
tions. What is more interesting is that, in the 
-trt-L correlation, the coefficients before 
the logarithms of -E pGa k and trt are almost exact 
opposites of each other. If this situation is input 
as a prior, the correlation is reduced to a 2D one, 
L oc (< cak /r RT )°- 842±0 - 064 , where E' peak and r RT 
denote the peak energy and minimum rise time 
in the GRB rest frame. We interpret/understand 
the result in the framework of the definition of the 
luminosity (energy released in units of time). Due 
to the complexity of GRBs, it is unlikely that we 
can accurately determine their luminosity through 
only one quantity, as can be seen from the large 
scatters of the 2D luminosity relations. And in 
principle, we should be able to calculate the lu- 
minosity from a characteristic energy scale and 
a characteristic timescale of GRBs. Our result 
suggests that -E pca k and trt are a pair of mea- 
surable quantities which are sufficiently correlated 
with the assumed characteristic energy scale and 
timescale and that there is a high degree of match 
between them (in the sense of the match between 
energy scale and timescale in the definition of the 
luminosity), so we may construct a better lumi- 
nosity relation with them. Since there are dif- 
ferent definitions of the luminosity, the index of 
0.842 in L oc (£ pcak /T RT ) a842±0 ' 064 may be inter- 
preted as the correction factor between the lumi- 
nosity attached to the pair of i? poa k and trt and 
the isotropic peak luminosity L. 

Our argument about the connection between 
the luminosity relations of GRBs and the defini- 
tion of the luminosity provides a clear direction 
for exploring tighter luminosity relations of GRBs 
in the future. It should also be easier to seek cor- 
relations between measurable quantities and the 
assumed characteristic energy scale and timescale 
related to the luminosity than to directly seek the 
relations between the measurable quantities and 
the luminosity. Our findings about the -E peak - 
trt—L correlation illustrate this approach. 

We thank the anonymous referee for helpful 
comments and suggestions. This work was sup- 
ported by the National Natural Science Founda- 
tion of China under grant no. 10973039, the Chi- 
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EW-W01, and the China Postdoctoral Science 
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